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Abstract. We have observed lES 1426+428 with INTEGRAL detecting it 
up to ~150 keV. The spectrum is hard, confirming that this source is an extreme 
BL Lac object, with a synchrotron component peaking, in a vF^, plot, at or above 
100 keV, resembhng the hard states of Mkn 501 and lES 2344+514. All these 
three sources are TeV emitters, with lES 1426+428 lying at a larger redshift 
[z — 0.129): for this source the absorption of high energy photons by the IR 
cosmic background is particularly relevant. The observed hard synchrotron tail 
helps the modeling of its spectral energy distribution, giving information on the 
expected intrinsic shape and flux in the TeV band. This in turn constrains the 
amount of the poorly known IR background. 



1. Introduction 

The spectral energy distribution (SED) of all blazars is characterized by two 
broad peaks, in a vFy plot, located in the IR/X-ray and in the MeV/TeV 
bands. The origin of the low energy peak is almost unanimously interpreted as 
synchrotron emission, while there is debate about the origin of the high energy 
peak. Coordinated and simultaneous variability (Maraschi et al. 1999) strongly 
suggests that most of the radiation of both components originates in one zone 
of the jet, and is due to the same particles. In the leptonic scenario the high 
energy peak is inverse Compton emission by the same electrons making syn- 
chrotron radiation at lower frequencies. These electrons can scatter the locally 
produced synchrotron photons (e.g., Maraschi et al. 1992) and any other impor- 
tant component of seed photons produced outside the jet, such as disk radiation 
(e.g., Dermer & Schlickeiser 1993; Celotti et al. 2007), photons coming from the 
broad line region (Sikora, Begelman & Rees 1994; Ghisellini & Madau 1996) or 
from the layer of the jet in a spine/layer model (Ghisellini et al. 2005). 

As suggested by Fossati et al. (1998) and Ghisellini et al. (1998) the 
blazar SEDs seem to form a sequence: in powerful blazars the high energy peak 
is dominant, and both peaks are at smaller frequencies. As the bolometric 
luminosity decreases, blazars becomes "bluer", both peaks shifting to higher 
frequencies (i.e. they become High Energy Peaked BL Lac, or HBL, in the 
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terminology of Padovani & Giommi 1995), with approximately equal power in 
the two peaks. At the low luminosity end of this sequence we then find BL Lacs 
which are TeV emitters. The synchrotron component of these BL Lacs peaks 
in the X-ray range, usually in the 0.1-1 keV range, but a few extreme BL Lacs 
have their synchrotron peak (at least occasionally), at or above '^lOO keV. 

lES 1426+428 was a candidate to belong to this extreme class: it was ob- 
served by BeppoSAX (Costamante et al. 2001) with the narrow field instruments 
(i.e. LEGS, MEGS between 0.1 and 10 keV) and the wide field PDS instrument 
(betwen 20 and 100 keV). The observed spectrum was hard, up to 100 keV, 
but there was some uncertainty concerning the PDS data, since this instrument 
could not resolve lES 1426+428 from the nearby source (another blazar) GB 
1428+4217. Figure [1] shows the BeppoSAX X-ray spectrum in which the PDS 
excess is clearly visible. The INTEGRAL observations have of course the re- 
quired angular resolution to disentangle the two sources, and have definitely 
confirmed the extreme nature of lES 1426+428. 

The X-ray spectrum of lES 1426+428 resembles the hard states of Mkn 501 
and 1ES2344+514, all belonging to the class of extreme BL Lacs, and all detected 
in the TeV regime (eg. Albert et al. 2007, Smith et al. 2006). The knowledge 
of the high energy synchrotron spectrum is crucial to construct a meaningful 
radiative model for the source, able to predict the intrinsic spectrum in the TeV 
band, produced (within the leptonic synchrotron self-Compton models) by the 
same electrons producing X-rays by the synchrotron process. This in turn is 
important to have information about the amount of absorption the TeV flux 
have suffered, due to the 7-7 —>■ process with the photons of the IR cosmic 
background. The relatively high redshift of lES 1426+428, z=0.129, makes it 
particularly interesting to study this effect. 
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Figure 1. The BeppoSAX spectra of 1ES1426+428. The spectral fit is 
derived from the low energy instruments, and extrapolated to the PDS range. 
The excess at high energies is evident. A similar excess is seen in the Bep- 
poSAX observation of GB1428+4217. See Costamante et al. 2001 for details 
of the fit. 
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2. INTEGRAL data 

A long look with INTEGRAL, 400 ksec, was given at this field in May 2006. 
For a description of the INTEGRAL mission refer to Winkler et al. (2003). The 
source was observed with the OMC (V filter), the JEM-X monitor (3-35 keV), 
the IBIS imager (0.02-10 MeV) and the SPI spectrometer (0.02-8 MeV). The 
source was not detected by SPI because it is below the detection limit. 

The blazar lES 1426+428 is detected, albeit at low significance, in both the 
JEM-X and the IBIS (ISGRI) instruments, while the quasar GB 1428+4217 is 
not detected in any instrument. We average the F-band light curve derived for 
lES 1426+428 from the OMC, deriving my=16.0. 

The SWIFT database (http://heasarc.nasa.gov/docs/swift/archive/) was 
searched to find observations closer in time than BeppoSAX to the INTEGRAL 
one. One of the longest observations, taken on May 19th, 2005, was used to 
derive an X-ray spectrum for lES 1426+428. 
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Figure 2. X-ray spectrum of lES 1426+428: from left to right the datasets 
are Swift /XRT (black), BeppoSAX/MECS (red), INTEGRAL/ JEM-X (blue) 
and INTEGRAL/ISGRI (green). 

Figure [2] shows the resulting X-ray spectrum, combining the Swift /XRT, 
BeppoSAX/MECS, INTEGRAL/ JEM-X and INTEGRAL/ISGRI datasets. 
JEM-X is the most uncertain: it needs to be better calibrated. In Table [U we 
report fluxes at different epochs, and a combined spectral shape from the same 
instruments as above. 

3. Discussion 

Figure [3] shows the SED of lES 1426+428, using most of the available datasets 
(see caption for details). The correction for the IR Extragalactic Background 
Light (IR EBL) is crucial for determining the correct intrinsic shape. We report 
two deabsorbed spectra, obtained using the "low IR" and "best" models of 
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Table 1. Table of fluxes and luminosities at different epochs with the dif- 



ferent instruments 



Instr. 


Date 


Flux (2-lOkGV) 

(10-11 cgs) 


Lx {2-lOkcV) 

(10^4 cgs) 


r [la range] 




XRT 


May 2005 


3.16 


13.9 


1.98 [1.97-2.00] 


1.32(392) 


MECS 


Feb 1999 


2.01 


8.8 






ISGRI 


May 2006 


1.58/1.65* 


7.0 







* the first number is the extrapolation to 2—10 kcV, the second is the measured flux in 20—100 kcV. 



Kneiske et al. (2004). The former case is compatible with the recent estimate 
of the background of Aharonian et al. (2006). 



— 8 I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — r 




Log u [Hz] 



Figure 3. The SED of lES 1426+428. Optical and ISGRI contemporaneous 
data as blue circles; MECS spectrum in red and other X-ray data (ASCA, 
ROSAT) in green; TeV observed points (CAT Djannati-Atai et al. 2002, 
HEGRA, Aharonian et al. 2003) in black, with two different de-absorption 
laws of the EBL (from Kneiske et al. (2004)): a milder one ("low IR") in 
blue, and a stronger one ("best") in magenta. 

The INTEGRAL spectrum shows that the peak of the synchrotron emission, 
even if in a relatively low observed state, is at energies higher than 100 keV. From 
ASCA data (Sambruna et al 1997) we know that the peak was in the 2-10 keV 
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range before 1999. All Swift observations, however, can be fitted (in the 2-10 
keV range) with a power law with F ~ 2 (Wolter et al. in prep.), indicating that 
the spectral shape did not change much recently, with flux variations of a factor 
2-3. Therefore we can reliably use these data to derive intrinsic quantities. 

3.1. Modeling the SED 
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Figure 4. The random Lorentz factor 7poak of the electrons emitting at 
the peaks of the SED as a function of the energy density (magnetic plus 
radiative, as seen in the comoving frame). The black dot corresponds to the 
lES 1426+428 modelled in this paper (see Figure[3]). Adapted from GhiseUini 
et al. 2002. 

We have applied a simple, one-zone, homogeneous synchrotron self-Compton 
model, described in Ghisellini, Celotti & Costamante (2002). The source is a 
cylinder of cross sectional radius R and width (in the comoving frame) AR' = 
R/T, where F is the bulk Lorentz factor. The magnetic field B is homogeneous 
and tangled. The particle distribution is the result of injection and cooling. We 
calculate the random Lorentz factor 7cooi at which the particles cool in one light 
crossing time AR' /c. Particles are injected between 71 and 72 with a power law 
distribution of slope s. In our case (7cooi > 7i) -^(7) oc 7"** between 71 and 
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7cooi! and A^(7) oc between 7cooi and 72. Below min(7cooi) 7i) we assume 

A^(7) oc 7""*^. We neglect any source of external photons. 

Based on the above assumptions, the SED modeling yields the results shown 
in Figure [3l For the adopted parameters, we find that in one light crossing time 
7cooi > 72, so that the particle distribution is A^(7) oc 7"^ up to 71 = 6000, where 
it steepens to A^(7) oc 7~2-86_ 'pj^jg electron slope makes a synchrotron spectrum 
F{u) oc v"^-^^, which is slightly increasing in a uFi, plot. As a consequence, 
the peak of the spectrum is made by the most energetic electrons, 7peak = 72 = 
3 X 10^. The magnetic field is B = 0.08 G, i? = 8 x 10^^ cm, F = 18, the viewing 
angle 9 = 1.8°, the Doppler factor 5 = 27.3. The power injected in relativistic 
electrons is L' = 3 x 10^^ erg s~^. 

These parameters are similar to the ones found for other low power, TeV 
emitting, BL Lac objects (see e.g. Katarzynski et al. 2006), even if 7peak is the 
largest ever found. 

In Figure U] we show 7peak as a function of the comoving energy density 
(magnetic plus radiative) as seen in the comoving frame, for a large sample of 
blazars. Note that the values found here for lES 1426+428 are along the well 
defined sequence, but reaching the extreme upper end in 7peak- 

4. Conclusions 

We observed lES 1426+428 with INTEGRAL: its spectrum is consistent with 
the BeppoSA'K observation, with a small variation (40%) in fiux. GB 1428+4217 
is not detected by INTEGRAL, in any band. 

A modeling of the SED in the standard SSC reference frame gives extreme 
parameters (the largest jpeak ever found), and the source confirms the physical 
interpretation of the blazar sequence. Our modeling prefers a moderate/low 
value of the cosmic background IR light. 
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